Vortices in a mesoscopic superconducting circular sector 
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In the present paper we develop an algorithm to solve the time dependent Ginzburg-Landau 
(TDGL) equations, by using the link variables technique, for circular geometries. In addition, 
we evaluate the Helmholtz and Gibbs free energy, the magnetization, and the number of vortices. 
This algorithm is applied to a circular sector. We evaluate the superconduting-normal magnetic 
field transition, the magnetization, and the superconducting density. Furthermore, we study the 
nucleation of giant and multi-vortex states for that geometry. 



I. INTRODUCTION 

The advances in the technologies of nanofabrication 
in the last few decades allowed intensive investigation 
efforts in nanostructured superconductors, both in the 
experimental and theoretical fronts. It is well known 
that, for very confined geometries, the superconducting- 
normal (SN) magnetic field transition is increased ex- 
traordinarily. It was experimentally observed that for an 
Al square very thin film with a size of few micrometers 
the upper critical field H C 2(T) can be increased up to 3.32 
with the inclusion of defects, 2.01 larger then the usual 
value of H C 2(T)^^ Furthermore, numerical simulations 
carried out in a circular wedge (see Ref. [3] and references 
therein), have shown that, by keeping the area of this 
geometry constant, the SN transition field is a uniform 
increasing function with decreasing the angular width 
and diverges as the angle goes to zero. 

Another important issue in confined geometries is the 
occurrence of giant vortices. The experimental obser- 
vation of giant vortex in mesoscopic superconductor is 
still a controversial issue. Through multi-small-tunnel- 
junction measurements in an Al thin disk film, Kanda et 
aZPhave argued that, as the vorticity increases, giant vor- 
tex configuration will occur. On the other hand, scanning 
SQUID microscopy on Nb thin film, both square and tri- 
angle, cannot guarant ee giant vortex configurations, at 
least for low vorticityJ^ Early numerical simulation of 
the present authors^ have shown the dynamic of the nu- 
cleation of giant and multi-vortex state before they set 
into an equilibrium configuration for a square geometry. 

It is well known that the phenomenology of super- 
conductivity can be described by the time dependent 
Ginzburg-Landau (TDGL) equationsP The present con- 
tribution uses the TDGL approach to address the issues 
above, namely, of the nucleation of vortices in confined 
geometries and the behavior of the transition field for a 
deformable geometry. For this, we have chosen a circular 
sector (see Fig. [T|), where we can arbitrarily change its 
shape. To our best knowledge, the discretization of the 
TDGL equations, by using the link variables technique, 
has been done only in rectangular coordinates! 10111 1 So, 
we will extend this algorithm to circular geometries by 



using polar coordinates. Our procedure makes possible to 
generalize the algorithm to any geometry. The key point 
in such problem is how to write the auxiliary fields ap- 
propriately according to the system of coordinates, mak- 
ing the development of the present algorithm a specific 
algorithm necessary. Otherwise, the purpose of general- 
ization will not be achieved. We anticipate that the two 
slopes of the present work will show that: (a) as we de- 
crease the area of the circular sector, the transition field 
may increase for large angles, but all the curves will col- 
lapse into the asymptotic behavior HjH^iT) = v3/@; 
and (b) only the confinement of vortices is not sufficient 
to obtain giant vortex state, but also the geometry is very 
important to favor the nucleation of such configurations. 
The geometry we have chosen does not allow us to assure 
clearly the nucleation of giant vortex. For a circular sec- 
tor of several angular widths, ranging from 45° to 180°, 
having the same area as a disk, a square and a trian- 
gle, we did not observe the occurrence of giant vortex as 
previous numerical simulations have predicted for those 
geometriesP In addition, we will show that the criterion 
used for nucleation of giant vortex may lead us to non 
conclusive pictures, at least for the geometry under the 
present investigation. 

The paper is outlined as follows. In Section[TT]we write 
the TDGL equations in a gauge invariant form by using 
the auxiliary field in polar coordinates. In Section III we 



develop the algorithm we use to solve the TDGL equa- 
tions: we define the mesh used to discretize the TDGL 
for a circular sector, the discrete variables which are eval- 
uated in the mesh, the boundary conditions and, finally, 
the important physical quantities which will be extracted 
from the numerical setup are determined. In Section |IV| 
we present and discuss the results of the numerical sim- 
ulations for certain parameters of a superconducting cir- 
cular sector. 



II. THE TDGL EQUATIONS 

The properties of the superconducting state are usually 
described by the complex order parameter ip, for which 
the absolute square value \tp\ 2 represents the superfluid 
density, and the vector potential A, which is related to 
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the local magnetic field as h = V x A. These quantities 
are determined by the TDGL equations, which in the 
non-dimensional version are given by 
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where T is the temperature in units of the critical temper- 
ature; lengths are in units of £(0), the coherence length 
at zero temperature, and fields in units of -ff C 2(0), the 
upper critical field at zero temperature; (3 is the ratio be- 
tween the relaxation times of the vector potential and the 
order parameter; k is the Ginzburg-Landau parameter 
which is material dependent; the operator D = — iV — A; 
Re indicates the real part of a complex variable and the 
overbar means the complex conjugation; (for more de- 
tails, see Reference I7I10I11[I . Here, we will neglect the 
z-dependence on the order parameter. This is valid ei- 
ther if the system is infinite along the z-direction or if it is 
a very thin film of thickness d <C 1. However, in the for- 
mer case, k 2 is replaced by an effective Ginzburg-Landau 
parameter k 2 ^ = n 2 /d (for instance, see References 
and fT3]) . The generalization to the a system of arbitrary 
thickness should not present any difficulty. 

It is convenient to introduce the auxiliary vector field 
U = (U p ,lAg) in polar coordinates, which is defined by 



U p (p,d) 



exp 



Ug(p,6) = exp f-i J 
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where (po,9q) is an arbitrary reference point. For the 
sake of brevity we omit the time dependence on the fields. 
Notice that 
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and that 
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Upon using these two last equations recursively, we 
obtain 
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As a consequence, we obtain for the kinetic term in the 
first TDGL equation 
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From equations Q, it can also be easily proved that 
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where Im indicates the imaginary part of a complex vari- 
able. 

Finally, by using equations ^ and ([7]), the TDGL 
equations of Q can be rewritten as 
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Disregarding the non-linear term, the first TDGL 
equation written as above resembles a diffusion equation, 
except by the fact that the Laplacian appears with dif- 
ferent weights. The weights depend locally on the com- 
ponents of the auxiliary field U. Written like in |8]) , the 
TDGL equations are gauge invariant, that is, they do not 
change their form under any transformation ip — > ipe lx , 
A — * A + V%. This is a very important point for any 
discretization procedure of the TDGL equations. Other- 
wise, we may obtain non-physical numerical solutions. 



III. NUMERICAL METHOD 

A. The Computational Mesh 

We will discretize the TDGL equations of ^ on a cir- 
cular sector as illustrated in Fig. [T] The mesh consists 
of N p x Ng cells with size (a p ,ag) in polar coordinates. 
The circular sector has internal radio r and external R; 
O is its angular width. Let (pi,9j) be a vertex point in 



the mesh, where pi+i = pi 



9j + i = 9j + ag, for all 



{1 < i < N p , 1 < j < Ng}; p\ — r and 9\ = 0; this partic- 
ular choice for the initial value of the angle docs not imply 
in lost of generality since the system is invariant under 
any rotation. The superconducting domain is compre- 
hended by fisc = {pi + a p /2 < p < p Np + a p /2, a e /2 < 
9 < 9n b + ag/2}. The superconducting region is sur- 
rounded by a thin normal metal layer of width a p /2 in 
the radial direction. Both regions are inside the domain 
O = {pi < p < pn p +i,0 < 9 < 9 Ng+ i}. We denote 
by dflsc the interface between the superconductor and 
the normal metal, and by d£l the normal metal-vacuum 
interface. 
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FIG. 1: The computational mesh in polar coordinates used 
for the evaluation of tpij (#, vertex point); h Zl ij and Li,j (Q, 
cell point); A p ,ij and U p> ij (■, link point); Ae,ij and Ue,i,j 
(□, link point). The superconducting domain is delimited by 
the dashed line dSlsc, and superconductor and the normal 
metal are surrounded by the solid line dQ.. Other details of 
the Figure are described in Section [HI A| 



B. The Discrete Variables 

Let us define the following discrete variables. 

• The vertex points 

Pl = r + (* - l)o p , 1 < * < N p + 1 , 

6j = (j-l)ofl, l<i<iV + l, (9) 

The points (p,:+i/ 2 = Pi + a p /2,0 i+1 / 2 = 0, t + ag/2) 
are called the cell points. The points (pi+i/2,#j) and 
(pi, Oj+1/2) are the link points in the radial and transver- 
sal directions respectively (see Fig. [lj . 

• The order parameter 



for all {1 < i < N p + 1 , 1 < j < Ng + 1}. 
• The vector potential 



(10) 



Ap t i,j = Ap(p i+1 / 2 ,0j) , A 8ti j — A 6 (pi,9j +1 / 2 ) , (11) 

for all {1 < i < N p , 1 < j < Ng + 1}, and {1 < i < 
N p + 1 ,1 < j < Ng}, respectively. 
• The link variables 

U p . hJ ^U p {p i ,Bj)Up{p i j rX ,d j ) = exp (-ia p A Piitj ) , 
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for all {1 < i < N p , 1 < j < N e + I}, and {1 < i < 
N p + 1 ,1 < j < Ng}, respectively. 
• The local magnetic field 

hz,i,j = h z (pi + i/2,0j + i/ 2 ) , 



(13) 



for all {1 < i < N p , 1 < j < N e }. 

In what follows, it will be important to define the fol- 
lowing discrete variable 



Li.j = exp — i 



A ■ dr 



exp [ ~i I h z pdpdO 
exp (-ia p p i+1/ 2a e h Zti j) , 



(14) 



for all {1 < i < N p , 1 < j < Ng}, where is the domain 
of a unit cell limited by a closed path d'S). The use of 
the Stoke's theorem and the midpoint rule for numeri- 
cal integration have been made. A simple inspection of 
equation (14) leads to 



Up^ijUe^+ijUpjj+iUg^^ 



(15) 



The Discretization of the TDGL 



Now, we are in a position to discretize the TDGL equa- 
tions. This can be done by using the central difference 
approximation for the derivatives which is second order 
accurate in (a p ,ag). A tedious, however straightforward 
calculation, leads us to the following discrete version of 
the TDGL equations of pl 
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From the numerical point of view, it is more convenient 
to evaluate the link variables rather than the vector po- 
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tential. From equations (12), we can easily verify that 
dA ;^ Uo.i.i dUni.j dA 



dt 
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In addition, from equation (14), we can write, accurate 
to second order in (a p ,ag) 



Im(l - Ljj) 
o-pPi+i/20-e 



(19) 



where Lij is given by equation (15). Upon introducing 



equations (18 1 and (19) into the second and third equa- 



tions of ( 16 1 we obtain the following recurrence relations 
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where 
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Finally, on using the one-step forward-difference En- 
ter scheme with time step At, we obtain the following 
recurrence relations 



^ id (t + At) = A, 3 (t) + At^ iW (t) 



U a ,ij(t + At) = Ua ti j(t) exp [ , 



(22) 



where a — (p, 9) . Notice that equations ( 22 1 were written 
in such a manner they guarantee the link variables are 
unimodular functions. The first recurrence relation run 
for all interior vertex points of f^sc, that is, {2 < i < 
N p ,2 < j < Ng}; the second ones for all link points in 
the interior of fl, that is, {1 < i < N p , 2 < j < Ng} for 
a = p, and {2 < i < N pi 1 < j < Ng} for a = 6. At 
the edge points of fi, the values of the discrete variables 
will be evaluated using the boundary conditions (sec next 
Section). 

There is a severe limitation on the choice of the time 
step At such that the recurrence relations converge. We 
have learned experimentally that the condition for sta- 
bility is assured by the following practical rule 



At < min 



5^ Pp 
4 ' 4k 2 



(23) 



where 



(24) 



Notice that the stability is controlled by the size of the 
smallest unit cell. The smaller the value of r, the more 
severe the restriction on the time step becomes. Perhaps 
in this case, it would be more convenient to use cither 



a semi or a full implicit scheme to solve equations ( 20 1 
which are usually unconditionally convergent. 



D. The Boundary Conditions 

Let n be a unit vector normal to the cT2gC inter- 
face and directed outward the domain f2sc- We will 
assume that the normal current density vanishes at the 
superconductor-normal metal interface, that is, • n = 
0. By using equations Q, it can be shown that the dis- 
crete implementation of this condition is as follows 



Up,N p ,j^N p , 3 
Ug t i t N a l]>i,Ng ■ 



(25) 



The first two equations run for all values of {2 < j < Ng}, 
and the second ones for all values of {2 < i < N p }. At the 
corner vertex points of the domain f2 are not necessary 
to run the recurrence relations (22 1. 



These last four equations update the values of the or- 
der parameter at any vertex point at the <9f2 interface. 
The values of the link variables at this interface will be 
updated by using the fact that the z-component of the 
magnetic field is continuous at the interface dflsci that 
is, h z> x,j = h z . Np , d = h Ztil i = h Z:l , Ng = H, which is 
the external applied magnetic field. Consequently, from 
equations (14 1 and (15 1, the link variables are updated 
according to 

L it j = exp (-ia p Pi + i/ 2 aeH) , (26) 
which runs for all edge points at the interface dfl. 

E. The Physical Quantities 

The topology of the superconducting state is usually il- 
lustrated by \ip\ 2 - This quantity can be determined from 
the outcome of the recurrence relations previously de- 
rived. Other important physical quantities used to de- 
scribe the vortex state are the Gibbs free energy, the 
magnetization, and the vorticity. In what follows we will 
derive an expression for each of these physical quantities. 
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• The condensation energy 
if c = (1-T) 2 JJ W 2 (^lV| 2 -l) Pdpdfl 
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• The field energy 

if £ = K 2 ff JJ h 2 zP dpde 
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The total Hclmholtz energy is then given by Jzf = Jz?k+ 
if c + iff. The Gibbs free energy can be obtained by a 
simple modification in the field energy. Instead of h 2 we 
would have (h z — H) 2 , or (Ji z ,i,j — H) 2 in the discrete 
version. Notice that the discrete TDGL equations could 
also be derived through the following equations 



(1-T) 



dt 



dA c 
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i as? 
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where srfi = a p pia s , £# p ,i = a p p l+1 / 2 a e , and £/ e ,i = 
a pPi a 6i which are the areas surrounded by the vertex and 
link points respectively. In order to derive the discrete 



TDGL equations by this means, it is essential to use the 
following relations 



dU 



dA c 



dU, 
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which can be easily shown from equations ( 12 ) 



The magnetization is 4nM = B — H, where B is the 
magnetic induction which is given by the spatial average 
of the local magnetic field. We have 



AttM - 



N p Ne 
i=l j=l 



H 



(32) 



where si is total area of the circular sector. 

The vorticity can be determined by integrating the 
phase ip in each unit cell of the mesh. We have, 



i,3 



1 

2tt 



N, 



Vcp-dr, N=^^„ t , j , 

>' i=lj = l 

where is a closed path with lower left and upper right 
corner at and (£ + l,j + 1), respectively. 

In all of our numerical simulations described in the 
next Section, we calculated N in order to make sure that 
the number of vortices agrees with what we see on the 
topological map of the order parameter. 



IV. RESULTS AND DISCUSSION 

The recurrence relations derived in the previous Sec- 
tion were implemented as follows. We started from the 
Meissner state, where ip = 1 and U p = Ug = 1 everywhere 
as the initial condition. Then we let the time evolves until 
the system achieves the stationary stated This is done 
by keeping the external applied magnetic field H con- 
stant. Next, we ramp up the applied field by an amount 
of AH. The stationary solution for H is then used as the 
initial state to determine the solution for H + AH, and 
so on. Usually we started from zero field and increased 
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FIG. 2: The magnetization curve as a function of the external applied magnetic field for four values of the angular width of the 
circular sector. Each jump in the magnetization indicates a phase transition. The corresponding configuration for each phase 
is indicated in Table [I] 



H until superconductivity is destroyed. As a criterion for 
termination of the simulation, we monitored the Gibbs 
free energy as a function of H . When the value of this 
quantity changes its sign, then the transition from the 
superconducting to the normal state sets in. 

The parameters used in our numerical simulations were 
K = 0.28, which is a typical value for AlP the other pa- 
rameters are d = 0.1, T = 0, and (5 = 1. The internal 
radio and the area of the circular sector were taken fixed 
for any value of angular width G. We used r = 1/ir and 
16-7T for the area, such that the external radio is given 
by R = ^32vr/e + r 2 . The reason for taking these pa- 
rameters as such is because it makes possible compari- 
son between our results and previous ones (see Ref. [8j 
and References therein) . The size of the mesh varied ac- 
cording to the value of O. As a criterion we have taken 
the length of the most external unit cell no larger then 
0.25 x 0.25. Since the order parameter varied most signifi- 
cantly over a distance £(T) (in real units), we are certain 
of not loosing this variation within this criterion. We 



ramp up the applied magnetic field, typically in steps of 
AiJ=10~ 3 . 

In Fig. [2] we present the magnetization versus external 
applied magnetic field curves for several values of the 
angular width. These pictures present a typical profile of 
a magnetization curve of a mesoscopic superconductor. 
It presents a series of discontinuities, in which each jump 
signals the entrance of more vortices into the sample. 
Notice that the lower critical field does not vary with the 
shape of the circular sector. From this, the immediate 
conclusion is that it depends only on the area but not 
on G. This is not an obvious result. In fact, in Ref. 8 
numerical simulations were performed in three different 
geometries: disk, square, and triangle and using the same 
parameters as in the present contribution. They find that 
the magnitude of the lower critical field is the same for the 
disk and the square, but slightly larger for the triangle, 
despite all geometries having the same area. Since we 
have the freedom to deform the circular sector, it should 
be expected similar behavior for the lower critical field, 
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FIG. 3: The nucleation field as function of the angular width 
of the circular sector. The solid line corresponds to equation 
( 34 1 taken from Ref . 3 and the open circles are the results 
found in the present simulation. 



especially for small width angles 0. However, the results 
shown in Fig. [2] do not seem to exhibit this feature. 

Another interesting feature present in the pictures of 
Fig. [2] is that the SN transition field H c3 {T) is approxi- 
mately the same for all angles greater than 90°. In addi- 
tion, HcsiT) is approximately the same for a disk and a 
square of equivalent area. However, for smaller values of 
O, this critical field becomes significantly larger. Indeed, 
in Ref. 3 H c3 (T) was calculated numerically for a wedge. 
The area they used for the wedge is 2.33 larger than the 
one we used here. They found that their results fit quite 
well into the following expression 



H c3 (T) 
H c2 (T) 



V3 

e 



1 + 0.148046' 



O.7460 2 
9 2 + 1.8794 I 



(34) 



In Fig. [3] we depict both the above expression and what 
we have found for the nucleation field 7J c3 (0) as a function 
of 6/77. As can be seen from that Figure, for large an- 
gles, the nucleation field is larger in the geometry we con- 
sider. This suggests that, had we diminished the area of 
the circular sector, this difference for large angles would 
have increased. Nonetheless, for small angles all curves 
should collapse into a single curve, which corresponds to 
the asymptotic behavior H c3 {T) / H c2 {T) = V3/Q deter- 
mined in Ref. [3] This suggests that the superconductor 
behaves as an unidimensional system as becomes small, 
no matter what the area is. 

We also investigated the topology of the order param- 
eter. Before going any further, let us establish the cri- 
terion we use to distinguish a single vortex from a giant 
vortex state. A giant vortex is nucleated as two or more 
vortices collapse into a single vortex in which all of them 
have, rigorously, a common core center. To describe an 
N vortex state we use the following nomenclature. We 



TABLE I: The sequence of vortex configurations for four dif- 
ferent angles of the circular sector. The nomenclature used is 
explained in the text. The configurations in brackets corre- 
sponds to what we obtain not using a logarithm scale. 
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1 sn° 


1 1^° 


90° 

•J\.) 


45° 


1 
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10 




1 ^ 

ID 
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4S 


4S 




4S 
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5S 


5S 


5S 
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6S 


6S 




6S 
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7S 


7S 


7S(5S1G 2 ) 
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8S 


8S(6S1G 2 ) 


8S(5S1G 3 ) 


8S(5S1G 3 ) 


9 


9S 


9S(4S1G 5 ) 


9S(5S1G 4 ) 


9S(4S1G 5 ) 


10 


10S(4S1G 6 ) 


10S(4S1G 6 ) 


10S(3S1G 7 ) 


10S(5S1G 5 ) 


11 


11S(5S1G 6 ) 


11S(2S1G 9 ) 


11S(3S1G 8 


11S(4S1G 7 ) 



denote by N S S, a multiple vortex configuration formed 
by N s single vortices. A single giant vortex of vorticity 
N g is denoted by IGn ■ For example, the ASIG2 state 
is formed by four single vortices and a double quantized 
giant vortex. 

In all geometries we have considered, usually it occurs 
transitions either from N to TV + 1 or from TV to N + 2 
vortices. In Fig. [4] and [5] we depict for = 180° 
and = 45° respectively, and two stationary states with 
different values of H . We have chosen transitions where 
we could have the formation of a giant vortex. As can be 
seen in these figures, we have the transitions 9S — > ASIG§ 
(0 = 180°) and 6S -> 551G 3 (0 = 45°) P However, if 
we look at the same pictures in a logarithm scale, we still 
see that the core centers of the vortices occupy different 
positions. So, within this criterion, we cannot affirm that 
a giant vortex has been nucleated. For higher vorticity, 
we have not observed any giant vortex either. A very dif- 
ferent scenario takes place in disks, squares, and triangles 
even using the same parameters as in the present workP 
Maybe, for smaller areas, giant vortex could be formed; 
we have not tested this possibility. All possible configu- 
rations are summarized in Table [I] up to N = 11. Notice 
that the vortices are always symmetrically distributed 
along the mediatrix. 



V. SUMMARY 

In summary, an algorithm has been developed for solv- 
ing the Ginzburg-Landau theory for circular geometries. 
This will probably make much easier to extend the tpU 
method for other geometries in addition to the circular 
an rectangular ones. Furthermore, we have applied the 
algorithm to circular sector and have found several con- 
figurations for the vortex state in this geometry. Also, the 
superconducting nucleation field has been evaluated. We 
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FIG. 4: (Color online) Two dimensional density contour plots of for an angular width of O = 180° with N = 9 and 
N = 10 vortices, from the top to the bottom (left column; the right column corresponds to the same pictures, but in logarithm 
scale). Both pictures correspond to the stationary states as the vortices enter the sample. Notice that, once the equilibrium 
configuration is achieved, the vortices are symmetric with respect to the vertical axis. This feature is always present in the 
other geometries. 

have presented some evidences that, as we diminishes the de Souza Silva for useful discussions. 

area of the supercondutor, the nucleation field increases. 

However, as the angular width goes to low values, this 

field exhibits an universal behavior, independently of the 

area. 



Acknowledgments 

The authors thank the Brazilian Agencies FAPESP 
and CNPq for financial support, and Dr. and Clecio C. 



1 G. R. Berdiyorov, B. J. Baelus, M. V. Milosevic, and F. 
M. Peeters PRB68, 174521 (2003). 

2 V. V. Moshchalkov et al. in Connectivity and Supercon- 
ductivity, edit by J. Berger and J. Rubinstein (Springer, 
Heidelberg, 2000. 

3 V. A. Schweigert and F. M. Peeters, Phys. Rev. B 60, 3084 
(1999). 

4 A. Kanda, B. J. Baelus, F. M. Peeters, K. Kadowaki, and 
Y. Ootuka, PRL93, 257002 (2004). 

5 T. Nishiro, S. Okayasu, J. Suzuki, and K. Kadowaki, Phys- 
ica C 412-414, 379 (2004). 

6 S. Okayasu, T. Nishio, Y. Hata, J. Suzuki, I. Kakeya, K. 
Kadowaki, and V. V. Moshchalkov, IEEE Trans. Appl. 



Supcond. , 678 (2005). 

7 E. Sardella, A. L. Malvezzi, P. N. Lisboa-Filho, and W. A. 
Ortiz, Phys. Rev. B 75, 014512 (2006). 

8 B. J. Baelus and F. M. Peeters, Phys. Rev. B 65, 104515 
(2002). 

9 A. Schmid Phsys. Kondens. Materie 5, 302 (1966). 

10 W. D. Gropp, H. G. Kaper, G. K. Leaf, D. M. Levine, M. 
Palumbo, and V. M. Vinokur, J. Comput. Phys. 123, 254 
(1996). 

11 G. C. Buscaglia, C. Bolech and A. Lopez, Connectivity and 
Superconductivity, J. Berger, and J. Rubinstein (Eds. ), 
Springer, 2000. 

12 J. Pearl, Appl. Phys. Lett. 5, 65 (1964). 



9 




FIG. 5: (Color online) The same as Fig. for 9 = 45° with N = 6 and N = 8. 



G. R. Berdiyorov, M. V. Milosevic, and F. M. Peeters, 
Physica C 437-438, 25 (2006). 

We use the follwoing criterion for the stationary state: if 
the highest difference ||V>(i)| — \ip (i+ At)|[, for any vertex 
point in the mesh, is smaller than a certain precision e, 
then we go over the next field. We have worked with a 
precision of e = 10~ 6 for 6 = 45° and e = 10 -5 for the 
other angular widths. 

We have used the hightest resolution as possible in order 
to detect an individuality in the vortex configurations. Be- 
yond some critical resolution, the pictures do not present 
an change. 



